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INTRODUCTION 


A  central  problem  in  developing  effective  smoke  obscurants 
for  space  applications  is  to  maximimze  volumetric  extinction  or 
absorption  at  the  wavelength  of  the  incident  electromagnetic 
radiation.  Yet  this  is  a  problem  without  a  general  solution. 
There  are  an  enormous  number  of  different  kinds  of  obscuring 
particles:  they  may  vary  in  composition,  shape,  and  size.  To 
choose  a  combination  of  particles  that  will  give  the  greatest 
obscuration  over  the  largest  possible  wavelength  interval  by 
doing  detailed  calculations  for  all  possible  particles--even  if 
the  necessary  analytical  tools  were  available--would  be  an 
endless  undertaking.  Nevertheless,  it  is  possible  to  make 
recommendations  based  on  very  general  considerations.  One  can 
establish  bounds  not  on  extinction  at  any  wavelength,  but  rather 
extinction  integrated  over  all  wavelengths,  by  invoking  sum 
rules. 

AIR  FCT-S  ^FTTCE  0?  SCIENTIFIC  RESEARCH  (AESC1 
5n?I  C£  OF  TRAN-CMI TTAL  TO  DTI  C 
This  technical  report  has  been  reviewed  and  is 
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Distribution  is  unlimited, 
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2.1  General  Considerations 


Sum  rules  follow  from  very  general  relations  called 
dispersion  or  Kramers-Kronig  relations.  Consider  two 
time-dependent  physical  quantities,  which  we  may  call  the  input 
and  the  output;  the  corresponding  Fourier  transforms,  which  are 
complex  functions  of  frequency  u>  ,  are  denoted  by  X^(w)  and 
X_(w).  If  the  relation  between  these  two  transforms  is  linear 


L  »->  U  »  i  »  «  u  t.l  |.l  »_*  ^Jl  fc-t  l.*1  |-»  Cf  W'  |>  t 


jW -. 


Xq(oj)  =  R(o»)Xi(aj), 


and  causal  (the  output  cannot  preceded  the  input  in  time),  then 


the  real  and  imaginary  parts  of  the  response  function  R(u>)  are 


connected  by  integral  relations  (subject  to  restrictions  on  the 


asymptotic  behavior  of  R) . 


Although  dispersion  relations  are  important  in  elementary 


particle  physics  (e.g.,  Goldberger,  1960)  and  some  areas  of 


optics,  particularly  optical  properties  of  solids  (e.g.,  Stern, 


1963),  they  have  been  largely  ignored  by  workers  in  the  field  of 


absorption  and  scattering  of  light  by  small  particles.  One 


exception  is  the  work  of  Purcell  (1969),  which  has  received 


little  attention.  Box  and  McKellar  (1978)  derived  an  expression 


for  integrated  extinction  similar  to  Purcell's  but  from  a 


completely  different  point  of  view.  Recently,  Purcell's  results 


were  rederived  and  extended  by  McKellar,  Box,  and  Bohren  (1982), 


and  an  alternative  treatment  was  given  by  Bohren  and  Huffman 


(1983,  p.  116).  This  treatment  is  outlined  below  because  it 


will  form  the  basis  for  further  investigation. 


Consider  a  single  particle  illuminated  by  a  plane. 


monochromatic  wave  of  frequency  u  and  wavenumber  k  propagating 


in  the  z  direction.  At  sufficiently  large  distances  r  from  the 


particle,  the  electric  field  E  scattered  by  it  is  transverse 


and  is  given  asymptotically  by 


E  % 
s 


ik  (r  -  z) 


-ikr 


kr  >>  1 


WWW 


(1) 


where  the  incident  field  (polarized  in  the  x  direction)  E-  is 
Ee„  and  E  =  E  exp(ikz  -  iut).  The  vector  scattering  amplitude  X 
in  the  forward  direction  determines  the  extinction  cross  section 
Cext  through  the  optical  theorem 

Cext  =  prRe{S(0)}'  (2) 

where  S(0)  is  the  projection  of  the  forward  scattering  amplitude 
onto  the  direction  of  polarization  of  the  incident  beam. 

Because  the  relation  (1)  is  linear  and  causual  (the 
scattered  wave  cannot  precede  the  incident  wave  in  time),  the 
real  and  imaginary  parts  of  S(0)  are  connected  by  integral 
relations.  However,  the  asymptotic  behavior  of  S(0)  is  such 
that  it  is  not  S(0)  itself  that  satisifies  these  relations  but 
rather  a  suitable  modification  of  it.  For  integral  relations  to 
exist,  S( 0 )  must  vanish  in  the  limit  of  infinite  frequency.  But 
it  follows  from  Fraunhofer  diffraction  theory  that  for 
sufficiently  high  frequencies 

S(0)  =  J2--  2G  ,  (3) 

4ttc 

where  c  is  the  speed  of  light  in  vacuo  and  G  is  the  geometrical 
cross  sectional  area  projected  onto  the  incident  beam  (for  a 
sphere  of  radius  a,  G  =  tt  a 2 )  -  The  function  F  formed  from  S 


vanishes  in  the  limit  of  infinite  frequency;  if  S(0)  is  analytic 
in  the  top  half  of  the  complex  frequency  plane  (this  is  required 
by  causality)  and  decreases  at  least  as  w2  at  low  frequencies, 
then  F  is  also  analytic,  hence  its  real  ( F f )  and  imaginary  (F'1) 
parts  are  connected  by  the  dispersion  relations 


F'  (w) 


SIF"  (fl) 

Q 2  -  oj2 


da  , 


(5) 


F  '  '  ( u) 


o 

2a) 

TT 


,00 


f  1  (a) 
a2  -  a)2 


da  , 


(6) 


o 

where  P  denotes  the  Cauchy  principal  value  of  the  integrals. 
If  we  now  substitute  (2)  in  (6),  we  obtain  the  dispersion 
relation  satisfied  by  the  extinction  cross  section 


P  — ext  (  -■  da  =  -2tt2c2  *  (7) 

J  a2  -  a)2  w 

0 

It  is  tempting  to  take  the  limit  of  (7)  as  the  frequency  goes  to 
zero  and  write 


lim  P 

U)-*-0 


c  .  (a) 

ext 


da  = 


a2  -  or 


c  (a) 

ext 


da  =  -2TT2c2lim 


2„2l Im{S (0) } 


a* 


U)->-0 


U) 


(8) 


0  0 

Indeed,  this  is  what  has  been  done  previously.  But  underlying 
(8)  are  some  implicit  assumptions  that  we  must  now  investigate. 

I 

i 


2.2  Requirements  on  the  behavior  of  C  .  at  low  frequencies 


Let  us  suppose  that  the  limit  of  the  right  side  of  (7)  as 


the  frequency  goes  to  zero  exists.  This  means  that  the  limit  of 
the  left  side  also  exists,  but  it  does  not  necessarily  follow 
that  this  limit  is 


C  (ft) 
ext 


dft 


ft  2 


(9) 


The  Cauchy  principal  value  is  also  obtained  as  a  limit,  and  the 
order  of  limits  cannot  always  be  interchanged.  Suppose,  for 
example,  that  Cext  were  a  constant.  Since  Cext  is  inherently 
positive,  it  follows  that  (9)  diverges  because  the  integral  of 
1/ft2  over  any  region  containing  the  origin  diverges.  Thus,  it 
is  clear  that  there  must  be  restrictions  on  the  low  frequency 
behavior  of  Cext  in  order  for  (8)  to  be  valid.  For  simplicity, 
let  us  assume  that  Cext  is  proportional  to  the  frequency  squared 
at  low  frequencies.  What  we  want  to  prove  is  that 


lim  P 

w->o 


ft2  -  U)z 


dft  » 
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Cext(n) 


dn  . 


(10) 


0  o 

From  the  definition  of  the  Cauchy  principal  value,  the  left  side 
of  (10)  is 


lim  lim 

5-»-o 


' 

C  (ft) 

ex  t 

dft  + 

ft2  -  u)2 
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C  „  (ft) 
ext 


dft 


ft 2  -  U)2 


(11) 
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C  „  (ft) 

ext 
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ft2  -  u2 


dft 


We  can  make  w  and  <5  as  small  as  we  wish  (but  not  zero)  subject 
to  the  requirement  that  u>>  5.  The  upper  limit  of  integration  6 


6 


can  also  be  made  as  small  as  we  wish  subject  to  the  requirement 
that  it  is  greater  than  w  +  6.  if  Cext  is  proportional  to  the 
frequency  squared  at  sufficiently  low  frequencies,  then  we  can 
substitute  the  low  frequency  expression  in  the  first  two 
integrals  in  (11).  These  integrals  can  now  be  evaluated,  and 
both  limits  taken.  The  result  is  (10).  Thus,  if  Cext  decreases 
as  the  frequency  squared  in  the  limit  of  zero  frequency,  (8)  is 
valid.  This  requirement  is  more  stringent  than  it  need  be.  All 
that  is  required  is  that  at  low  frequencies  Cgxt  decrease  as 
frequency  to  the  power  1  +  a  ,  where  o  is  any  positive  number. 

At  sufficiently  low  frequencies,  the  extinction  cross 
section  of  an  ellipsoid  composed  of  either  a  simple  free 
electron  metal  or  a  simple  insulator  is  proportional  to 
frequency  squared  (see,  e.g.,  Bohren  and  Huffman,  1983,  pp.  345 
and  348).  Thus,  the  low  frequency  behavior  of  Cext  for  real 
particles  is  such  that  (8)  is  valid. 

It  also  follows  from  the  preceding  analysis  that  there  do 
not  exist  sum  rules  for  the  extinction  cross  section  with  any 
form  other  than  (9).  For  example,  suppose  we  inquire  into  the 
existence  of  a  sum  rule  for  the  quantity  wCext*  The  low 
frequency  limit  of  this  quantity  is  such  that  we  may  interchange 
the  order  of  limits;  that  is. 


,00 

lim  p  — 

U)->-0 


ftC  (fl) 

ext 


n 2  -  oj 2 


dn 


Cext(f]) 

Q 


d  n. 


0  0 

Now,  however,  cext  does  not  have  the  proper  behaviour  at  high 
frequencies;  that  is,  <uF  does  not  vanish  in  the  limit  of 


7 


\ 

t 

t 

i 

i 

( 


infinite  frequency,  which  is  required  for  the  validity  of  the 

integral  relations  (5)  and  (6).  The  same  argument  holds  for  any 

quantity  of  the  form  '-AnCext,  where  n  is  any  number  other  than 

0.  Thus,  we  conclude  that  only  one  sum  rule  for  C  .  exists. 

ext 


2.3  Sum  rules  when  magnetic  dipole  terms  are  included 


To  proceed  further  we  need  to  determine  the  limit  on  the 
right  side  of  (8),  that  is,  we  need  to  know  the  low-frequency 
limit  of  the  forward  scattering  amplitude.  For  particles 
composed  of  an  electrical  insulator,  this  amplitude  will  be  that 
given  by  electrostatics,  i.e.,  the  electric  dipole  term  in  a 
multipole  expansion  of  the  electronic  ^netic  field  sec.  ctered  by 
the  particle.  But  for  metallic  particles,  it  is  not  clear  if 
this  choice  is  correct:  at  far  infrared  frequencies,  the 
magnetic  dipole  term  can  be  equal  to  or  even  greater  than  the 
electric  dipole  term.  For  this  reason,  the  role  of  the  magnetic 
dipole  term  in  sum  rules  was  investigated. 

In  order  to  apply  Rayleigh  theory  (i.e.,  electrostatics)  to 
scattering  by  small  particles,  two  conditions  must  be  satisfied: 
the  size  parameter  x  =  2iTa/X  ,  where  a  is  a  characteristic 
linear  dimension  of  the  particle,  must  be  small  compared  with 
unity  as  must  |m|x,  where  m  is  the  particle’s  complex  refractive 
index.  Let  us  consider  a  simple  f ree-electron  metal,  the 
dielectric  function  of  which  is  given  by  the  Drude  formula 

CO2 

e  =  e'  +  ie"  =  1  -  - 2 -  ,  (12) 

w2  +  iyw 


•-V  ' 


o 


where  is  the  plasma  frequency  (typically  in  the  ultraviolet) 
and  Y  is  a  damping  factor,  which  also  has  the  dimensions  of 
frequency.  The  refractive  index  is  the  square  root  of  e  .  The 
size  parameter  x  can  be  written  as 

x  =  —  ,  (13) 

from  which  it  is  clear  that  x  vanishes  in  the  limit  of  zero 
frequency.  But  since  |m|  does  not  vanish  in  this  limit  for  a 
material  described  by  (12),  it  is  not  obvious  that  |™|xdoes. 

The  modulus  of  the  refractive  index  is 

|  m  |  =  (t'z  +  e"2)*  .  (14) 

For  frequencies  oj  much  less  than  y  ,  and  subject  to  the 
requirement  that  Y<<  ,  we  can  combine  (12) -(14)  to  obtain 

a)  a  _ 

!m|x  =  -  /oj/y  o)  <<  y  <<  o)  (15) 

c  p 

Thus,  provided  that y  is  not  identically  zero,  |m|x  also  vanishes 

in  the  limit  of  zero  frequency.  This  condition  on  y  is 

equivalent  to  requiring  the  dc  conductivity  o>2/y  to  be  finite. 

p 

To  proceed  further,  we  must  consider  a  specific  particle. 

The  forward  scattering  amplitude  for  a  sphere  is  an  infinite 

series  in  the  scattering  coefficients  a  and  b  : 

n  n 


S(0) 


^(2n  +  1 )  ( a_  +  b  ) 


(16) 


For  sufficiently  small  x  and  |m|x,  the  electric  dipole  term  a^ 
and  the  magnetic  dipole  term  dominate  the  series  (16)  and  are 
given  approximately  by 


a 


1 


•  2x3  e  -  1 
3  e  +  2 


(17) 


If  y  is  not  identically  zero,  then  for  a  material  with 
dielectric  function  (12) 

lim  b, /a.  =0  (18) 

On  physical  grounds,  we  expect  (18)  to  be  valid  for  particles  of 
arbitrary  shape.  Thus,  we  have  established  that  the  sum  rule 
(8)  is  valid  even  for  metallic  particles  and  that  the  scattering 
amplitude  in  the  zero  frequency  limit  corresponds  to  electric 
dipole  radiation,  provided  only  that  the  dc  conductivity  is 
finite . 

Before  proceeding,  we  need  to  consider  the  consequences  of 
mean  free  path  limitations  for  sum  rules. 

2.4  Limitation  of  the  mean  free  path 

It  will  again  be  useful  to  consider  a  specific  particle,  a 


sphere.  We  also  define  the  volumetric  attenuation  coefficient 
^  v  as  the  extinction  cross  section  per  unit  particle  volume.  It 
follows  from  (8),  (16),  and  (17)  that  the  volumetric  attenuation 


coefficient  is  independent  of  particle  size  and  depends  only  on 
the  low  frequency  limit  of  the  dielectric  function 


r 

j  ot  dA  =  3  rr 2  lim  Re 

!  v 
J 


(  19) 


We  have  also  transformed  the  variable  of  integration  from 
frequency  to  wavelength.  For  a  material  with  dielectric 
function  (12),  the  term  on  the  right  side  of  (19)  is 


Re 


e  (w) 
,£  (w) 


-  n 

+  2 


-a) 2  (  3oj  2  -  oj 2 ) 

- E - E -  . 

(  3u)2  -  uj2  )  2  +  9y  2 u) 2 

P 


(20) 


The  limit  of  (20)  as  the  frequency  goes  to  zero  is  unity, 
regardless  of  the  value  of  y,  and  in  this  instance 


co 

a  dA  =  3tt 2  .  (21) 

v 

o 

The  damping  term  y  arises  from  scattering  of  electrons  by 
the  lattice,  electron-phonon  scattering.  If  it  were  not  for 
electron-phonon  scattering,  the  dc  conductivity  of  a  metal  would 
be  infinite.  In  a  sufficiently  small  particle,  scattering  of 
electrons  by  its  boundary  may  be  more  important  than  scattering 
by  the  lattice.  That  is,  the  mean  free  path  for  scattering  by 
the  boundary  may  be  less  than  that  for  scattering  by  a  phonon. 
This  was  first  pointed  out  by  Doyle  (1958),  and  subsequently 
investigated  by  Doremus  (1964),  Kreibig  and  von  Fragstein 
(1969),  Kreibig  (1974),  and  Granqvist  and  Hunderi  (1977). 

Kreibig  (1974)  interpreted  absorption  by  small  silver  particles 


by  modifying  the  damping  constant  in  the  Drude  expression  (12). 
This  damping  constant,  which  is  the  inverse  of  the  collision 
time  for  conduction  electrons,  is  increased  because  of 
additional  collisions  with  the  boundary  of  the  particle.  Under 
the  assumption  that  the  electrons  are  diffusely  reflected  at  its 
boundary,  can  be  written 


Y 


Ybulk 


+  Vp/L, 


(22) 


where  Ybulk  is  the  bulk  metal  damping  constant,  vp  is  the 
electron  velocity  at  the  Fermi  surface,  and  L  is  the  effective 
mean  free  path  for  collisions  with  the  boundary.  Kreibig  took  L 
to  be  4a/ 3  for  a  sphere  of  radius  a.  We  have  already  shown, 
however,  that  (21)  is  independent  of  Y  .  Thus,  we  conclude  that 
the  limitation  of  the  electron  mean  free  path  by  the  boundaries 
of  a  small  particle  do  not  affect  its  integrated  extinction. 
Detailed  calculations  for  spheres  (see  Bohren  and  Huffman,  1983, 
p.  338)  are  in  accord  with  this  conclusion.  When  particles  are 
sufficiently  small,  the  effect  of  their  boundaries  is  to 
decrease  peak  extinction  while  simultaneously  broadening  it  in 
such  a  way  that  integrated  extinction  remains  constant. 


3.  SUM  RULES  FOR  PARTICLES  OF  VARIOUS  SHAPES 


3.1  A  general  result  and  its  interpretation 


The  scattering  amplitude  S(0)  for  any  particle  that  is 


L2 


sufficiently  small  compared  with  the  wavlength  of  the  light 
illuminating  it  is  proportional  to  its  volume.  Thus,  it  follows 
from  (8)  that  integrated  extinction  per  unit  volume  of  any 
particle  is  independent  of  its  size:  it  depends  only  on  its 
composition.  When  a  particle  is  illuminated  by  an 
electromagnetic  wave,  various  electromagnetic  modes  may  be 
excited.  The  larger  the  particle,  the  greater  tha  number  of 
such  modes.  This  follows  from  simple  physical  reasoning, 
although  it  is  also  evident  from  exact  solutions  to  scattering 
problems.  The  smaller  the  particle,  the  fewer  the  modes  that 
can  be  excited.  Yet  integrated  extinction  is  independent  of 
particle  size.  That  is,  a  fixed  amount  of  total  extinction  must 
be  shared  by  modes,  the  number  of  which  increases  with 
increasing  particle  size.  Thus,  the  larger  the  particle,  the 
more  a  fixed  amount  of  total  extinction  has  to  be  shared  with  an 
ever  larger  number  of  modes.  Hence  the  extinction  associated 
with  each  resonant  mode  must  decrease  with  increasing  size. 

This  is  turn  implies  that  the  smallest  particles  have  the 
greatest  extinction.  This  is  a  very  general  result.  It  can  be 
verified  in  specific  instances.  For  example,  the  greatest 
possible  extinction  per  unit  volume  of  a  sphere  is  obtained  for 
spheres  much  smaller  than  the  wavelength.  This  high  volumetric 
extinction,  however,  is  concentrated  in  a  narrow  spectral  band, 
as  implied  by  the  sum  rule. 


3.2  Ellipsoidal  particles 


Equation  (8)  is  completely  general.  To  apply  it  to  a 
definite  particle  requires  the  scattering  amplitude  in  the 
forward  direction  in  the  limit  of  zero  frequency.  Consider  a 
homogeneous  ellipsoidal  particle  of  volume  v  illuminated  by 
light  polarized  along  one  of  its  principal  axes.  At 
sufficiently  low  frequencies,  this  particle  is  equivalent  to  a 
dipole  with  (induced)  moment 


P 

w 


£  V 
m 


£  -  1 _ 

1  +  L(£  -  1) 


(23) 


where  EQis  the  incident  field,  L  is  a  geometrical  factor  lying 
between  0  and  1,  and  em  is  the  permittivity  of  the  surrounding 
medium  (assumed  to  be  nonabsorbing).  In  general,  there  is  a 
distinct  geometrical  factor  Lj  for  each  of  the  three  principal 
axes  of  the  ellipsoid,  the  sum  of  which  is  unity.  In  general, 
the  incident  light  is  not  polarized  along  one  of  the  principal 
axes,  so  we  have  to  resolve  the  incident  electric  field  into  its 
components  along  the  principal  axes  of  the  elliposoid.  We 
denote  these  coordinate  axes,  which  are  fixed  relative  to  the 
particle,  by  xyz.  Axes  relative  to  which  the  incident  light  is 
specfied  (laboratory  coordinate  system)  are  denoted  by  x'y'z'. 

If  the  incident  light  is  polarized  along  the  x'  axis,  the 
scattering  amplitude  corresponding  to  (23)  is 
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where  e  is  a  unit  vector  in  the  scattering  direction.  The 
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projection  of  (24)  onto  the  x1  axis  in  the  forward  direction  is 


S(0) 
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where 


*11  aiail  +  a2a21  +  a3a31' 


a  =  v 


e  ~  1 


1  +  L  (e  -  1) 


(25) 


The  quantities  a^  etc.  are  direction  cosines  between  the  primed 
and  unprimed  axes,  and  they  satisfy 
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If  we  combine  substitute  (25)  in  (18)  we  obtain  the  sum  rule 

_S5£— -  dSi  -  lim  im 
Ci2  2c 

0 

As  a  check  on  (26),  we  note  that  it  reduces  to  (19)  when  all  the 
geometrical  factors  are  equal  (i.e.,  the  ellipsoid  is  a  sphere) 
and  the  variable  of  integration  is  changed  from  frequency  to 
wavelength. 

A  less  general,  but  more  useful,  form  of  (26)  is  obtained 
by  averaging  over  all  orientations  of  the  ellipsoid.  Unless 
special  care  is  taken,  a  cloud  of  ellipsoidal  particles  is 
likely  to  be  randomly  oriented,  in  which  instance 
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(26) 


where  the  brackets  indicate  an  average  over  all  directions  with 


equal  probability  per  unit  solid  angle.  The  corresponding 
volumetric  attenuation  coefficient  satisfies 


,00 


a  dX  =  7T2  lim  Re 
v 

u>-*-o 


e  -  1 


j=l  3  +  3Lj (e  -  1) 


(27) 


where  a  =  <cext>/v*  Again,  this  reduces  to  (19)  when  all  the 


geometrical  factors  are  equal. 

Each  term  in  the  sum  on  the  right  side  of  (27)  has  the  form 


( e '  -  1)  (3  +  3L  e'  -  3L.)  +  3L.e">2 

lim  - ^ - 2 - 1 -  (28) 

gj+o  (3  +  3L  e'  -  3Lj  )  2  +  9L*e"2 

Let  us  consider  the  low  frequency  limit  of  (28)  for  two  classes 
of  materials,  metals  and  polar  insulators.  The  Drude  dielectric 
function  (28)  applies  to  an  ideal  free-electron  metal.  The 
imaginary  part  of  this  dielectric  function  goes  to  infinity  in 
the  limit  of  zero  frequency.  Hence,  for  an  ideal  free-electron 
metal  ( 27 )  becomes 


,co 


a  dX 
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(29) 
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Now  consider  an  insulator,  one  for  which  the  imaginary  part  of 
the  dielectric  function  goes  to  zero  in  the  low-frequency  limit. 
For  such  a  material  (28)  is 
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lim  _ -1^ - 

jj-^o  3  +  3L .  £  —  3L . 

1  J 


Thus,  if 


L.  >>  lim  1/ ( e  •  -  1)  (j  =  1,2,3), 

-1  jJ->-0 


integrated  volumetric  extinction  is  still  given  (approximately) 
by  (29). 

All  the  geometrical  factors  L ^  for  a  sphere  are  1/3.  The 
low-frequency  limit  of  the  dielectric  function  of  water  is  more 
than  80.  Thus,  we  obtain  the  unexpected  result  that  integrated 
extinction  by  a  water  droplet  is  (approximately)  that  same  as 
that  by  an  equal-volume  metallic  sphere.  The  extinction  spectra 
of  these  two  particles  are  markedly  different,  and  yet  the 
integrals  of  these  spectra  are  the  same.  This  similarity  is  not 
maintained  for  all  shapes.  As  the  particle  shape  deviates  more 
from  sphericity,  at  least  one  of  the  geometrical  factors  must 
decrease  to  the  point  where  (31)  is  no  longer  satisfied. 
Nevertheless,  for  a  wide  range  of  shapes  integrated  (volumetric) 
extinction  by  metallic  particles  is  approxmately  equal  to  that 
of  insulating  particles  made  of  a  substance  with  a  high 
dielectric  at  low  frequencies  function  (e.g.,  water  in  both 
liquid  and  solid  phases). 


4.  SOME  THOUGHT  ON  OPTIMIZING  ATTENUATION 


The  fraction  T  of  electromagnetic  radiation  transmitted  by 


a  cloud  of  particles  is  a  product 


T  =  II  exp  (-a  .  f  .  z) 
j  VJ  J 


where  a  is  the  volumetric  attenuation  coefficient  of  the  jth 


kind  of  particle,  f ^  is  the  volume  fraction  in  the  cloud  (the 


subscript  j  may  denote  composition,  size,  shape,  or  even 


orientation),  and  z  is  the  cloud  thickness.  Equation  (32)  is  a 


lower  limit  on  the  cloud  transmission;  T  can  be  greater  because 


of  multiple  scattering.  The  quantity  f^z  is  the  total  volume  of 


particles  of  the  jth  kind  per  unit  area  of  the  cloud. 


Transmisison  can  therefore  be  made  arbitrarily  low  by  increasing 


the  mass  of  particles,  which  is  effective  but  not  efficient.  To 


minimize  the  amount  of  material,  one  must  maximize  . ,  and  not 


merely  over  a  narrow  range  of  wavelengths.  We  can  rewrite  (32) 


-InT/fz  =  E 


(33) 


where 


j  fJ  '  f3  '  V£  '  \  f3 


So  we  want  to  maximize  the  sum  on  the  right  side  of  (33).  I 


argued  in  the  Introduction,  that  to  do  this  by  detailed 


calculations  would  be  an  endless  task,  even  if  all  the 


analytical  tools  were  at  hand  (which  they  are  most  definitely 


not).  Nevertheless,  one  can  establish  some  bounds. 

First,  we  can  establish  some  practical  lower  bounds.  The 
extinction  cross  section  of  any  particle  larger  than  the 
wavelength  of  the  light  illuminating  it  is  approximately  2Q, 
regardless  of  its  composition,  where  G  is  the  geometrical  cross 
sectional  area  of  the  particle  projected  onto  the  incident  beam. 
G  is  proportional  to  some  characteristic  length  d  of  the 
particle  (for  a  sphere,  its  diameter).  The  volume  of  the 
particle  is  proportional  to  the  cube  of  d.  Thus,  for  such  a 
particle 

°v  =  Cext/V  =  K/d'  (34) 

where  K  is  a  constant  depending  on  particle  shape  (K  =  3  for  a 

sphere).  Let  us  take  d  to  be  of  order  1  urn,  which  gives  av  of 
4  -1  .  . 

order  10  cm  ;  if  d  is  ten  times  larger,  then  «v  is 

correspondingly  less.  Thus,  d  ^  1  Mm  is  more  or  less  the 

smallest  size  for  which  (34)  is  a  good  approximation  for  visible 

and  near-visible  wavelengths.  It  therefore  seems  that  a 

4  -1 

volumetric  attenuation  coefficient  of  order  10  cm  represents 
a  kind  of  lower  bound:  it  is  the  value  that  can  be  obtained 
without  making  much  of  an  effort.  What  is  the  upper  bound? 

According  to  the  sum  rule  (27)  integrated  volumetric 
attenuation  is  independent  of  particle  size.  I  argued  in  Sec. 
3.1  that  when  a  particle  is  illuminated  by  an  electromagnetic 
wave,  various  electromagnetic  modes  may  be  excited.  A  specific 
example  of  this  is  provided  by  a  sphere.  The  extinction  cross 


section  of  a  sphere  is  given  by  the  infinite  series 


Cext  =  ( 2^/k2 )Re{ £  (2n+l)(an  +  bn) }  (35) 

Each  coefficient  corresponds  to  an  electromagnetic  mode  of  the 
sphere;  an  correspond  to  modes  of  electric  type,  bn  to  modes  of 
magnetic  type.  The  larger  the  sphere,  the  more  terms  in  (35) 
are  needed  for  satisfactory  convergence;  that  is,  the  larger  the 
sphere,  the  more  modes  are  excited  by  the  incident  wave. 
According  to  the  sum  rule  (27),  the  integrated  (volumetric) 
extinction  is  independent  of  particle  size.  Thus,  the  larger 
the  particle,  the  more  that  a  fixed  amount  of  total  extinction 
has  to  be  shared  by  more  modes.  Conversely,  the  smaller  the 
particle,  the  fewer  the  modes,  hence  the  more  extinction 
associated  with  each  one  of  them.  To  make  my  point,  I  have 
invoked  a  sphere,  but  my  conclusion  is  independent  of  particle 
shape . 

We  therefore  must  look  to  small  (compared  with  the 
wavelength)  particles  for  the  greatest  possible  extinction  per 
unit  volume.  It  follows  from  (27)  and  (29)  that  for  a  given 
shape,  metallic  particles  will  have  the  greatest  integrated 
extinction.  Moreover,  it  is  the  extreme  shapes  that  gives  (29) 
its  greatest  value.  For  a  (perfectly  sharp)  needle  or  a 
(perfectly  flat)  disc,  =  0,  hence  for  these  ideal  particles 
integrated  extinction  is  infinite. 

To  estimate  the  greatest  possible  volumetric  extinction,  I 
did  some  simple  calculations  for  very  small  (compared  with  the 


2 


K 


I 

i 

wavelengths  of  ultraviolet  light)  spheres  and  randomly  oriented 

discs  and  needles  of  aluminum.  Aluminum  is  almost  an  ideal  ! 

I 

f ree-electron  metal;  its  dielectric  function  is  given  to  good  < 

i 

approximation  by  (12).  I  took  the  plasma  frequency  (in  electron  j 

volts)  to  be  15  eV;  the  damping  factor  was  0.6  eV.  The  results  1 

are  shown  in  Fig.  1.  Note  that  the  sphere  gives  the  greatest 
peak  extinction.  This  is  not  surprising  in  light  of  the  sum 
rules.  Fewer  modes  are  excited  in  a  sphere  than  in  an  < 

asymmetric  particle  such  as  a  needle  or  a  disk.  Peak  extinction 
is  about  107  cm-'*'.  I  claim  that  this  is  a  practical  upper 

i 

bound.  Moreover,  this  bound  cannot  be  obtained  over  a  wide 
range  of  frequencies.  A  practical  upper  bound  over  at  least  one 
decade  from  the  ultraviolet  to  the  near  infrared  is  106  cm-1, 
and  this  would  have  to  be  obtained  with  a  mixture  of  particles 
having  different  shapes. 

I  argued  at  the  beginning  of  this  section  that  volumetric 

4-1  I 

extinction  of  10  cm  could  be  obtained  with  almost  no  thought. 

One  could  almost  select  particles  at  random.  With  a  bit  of 

thought  it  is  possible  to  increase  this  number  by  a  factor  of 

100.  A  further  factor  of  ten  increase  might  be  possible,  but 

this  would  likely  be  obtained  at  the  expense  of  high  extinction 

in  one  spectral  region  being  paid  for  by  much  lower  extinction 


in  another. 


RELATIVE  WAVELENGTH 
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